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Abstract 

We investigate the emergent universe scenario in the presence of interacting fluids. 
The non-linear equation of state (EoS) considered in the general theory of relativity 
for obtaining emergent universe is effectively a cosmological model with a composition 
of three fluids. In this paper we consider two models to realize viable cosmological 
scenarios, viz., (i) a two-fluid model with interaction of a pressureless fluid with the 
fluid having the non-linear EoS needed for the emergent universe, and (ii) a three- 
fluid model with interaction among the three fluids which originate from the EoS of 
the emergent universe. It is found that realistic cosmological models in accordance 
with observations are not ruled out for both the above cases. We further show that 
the generalized second law of thermodynamics is found to hold good in the emergent 
universe with interacting fluids. 

PACS number(s) : 04.20.Jb, 98.80.Cq, 98.80.-k 


‘Electronic mail : bcpaul@iucaa.ernet.in 
^Electronic mail : archan@bose.res.in 


1 



1 Introduction 


The astronomical and cosmological observations predict that we live in an expanding Uni¬ 
verse. After the discovery of the cosmic microwave background radiation (CMBR) [1], the 
big-bang cosmology has become the standard model for cosmology which has a beginning 
of the Universe at some finite past. It is found that the big-bang cosmology does face some 
problems in addressing issues of the observed universe both in the early and late universe. 
A number of problems, namely, horizon problem, flatness problem, singularity problem etc. 
crop up when one probes the early universe in the framework of the big-bang model. These 
problems however can be resolved by evoking a phase of inflation la El n El El El El [9] at 
a very early epoch. Furthermore, it is known that the large scale structure formation of 
the universe can be addressed in this scenario successfully. On the other hand, very recent 
observations predict that our universe is passing through a phase of acceleration uniEiiEa. 
This phase of acceleration is believed to be a late time phase of the universe and may be 
accommodated in the standard model by adding a positive cosmological constant in the 
Einstein’s field equation. 

In spite of its overwhelming success, modern big-bang cosmology still has some unre¬ 
solved issues. The physics of the inflation and the introduction of a small cosmological 
constant for late time acceleration are not completely understood naEn in details. More¬ 
over, various competing models exist that are as yet not fully distinguished empirically from 
the currently available observational data. This is why there is enough motivation to search 
for an alternative cosmological model. In this context, Ellis and Maartens |15] considered 
the possibility of a cosmological model [l6] in which there is no big-bang singularity, no 
beginning of time, and the universe effectively get rid of a quantum regime for space-time 
by staying large at all times. The universe started out in the infinite past in an almost 
static Einstein universe, and subsequently, it entered in an expanding phase slowly, even¬ 
tually evolving into a hot big-bang era. Later Ellis, Murugan and Tsagas m constructed 
an emergent universe scenario for a closed universe with a minimally coupled scalar field 
4>, which has a special form for the interaction potential V{(j)). It was pointed out later 
[18j that the potential is similar to what one obtains from a modified gravitational action 
with a polynomial Lagrangian L = R-\- aB? after a suitable conformal transformation and 
identifying the field as cp = —\/3 ln(l -|- 2aR) with a negative a. 

The EU scenario merits attention, as it promises to solve several conceptual as well as 
technical issues of the big-bang model. A notable direction is regarding the cosmological 
constant problem |19] . Mukherjee et al. [20] obtained an emergent universe in the framework 
of general theory of relativity in a flat universe with a nonlinear equation of state of the 
form 

p = Ap-By/^ ( 1 ) 

where A and B are arbitrary constants. Such potentials have been also studied in the context 
of modified Chaplygin gas models m- A spatially flat universe is most likely as predicted 
from recent cosmological and astronomical observations. Such an emergent universe scenario 
can be realized also in a modified theory of gravitation m by including a Gauss-Bonnet 
term in the presence of dilaton coupling [23] , Brane world gravity [291 EHl 126] , Brans-Dicke 
theory m, the non-linear sigma model [28], Chiral cosmological fields in Einstein Gauss- 
Bonnet gravity [29], dark sector fields in a chiral cosmological model [30] and exact global 
phantonical solution m- Emergent universe accommodates a late time de-Sitter phase and 
thus it naturally leads to the late time acceleration of the universe, as well. Such a scenario 
is promising from the perspective of offering unified early as well as late time dynamics 
of the universe [32]. Note, however, that the focal point of unification in such emergent 
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universe models lies in the choice of the equation of state for the polytropic fluid, while 
several other models of unification rely more on the scalar field dynamics through choice of 
field potentials [33l EH ES] • 

The emergent model proposed by Mukherjee et al. m in which a polytropic equation of 
state (henceforth, EoS) is used, gives rise to a universe with a composition of three different 
types of fluids determined by the parameters A and B. In the original emergent universe 
model proposed by Mukherjee et al. [20], it was assumed non-interacting fluids and each of 
the three types of fluids identihed satisfy conservation equations separately. Recently using 
the observational prediction of WMAP7 [36] and Planck2013 m the permitted range of 
values of the parameters A and B are determined [38l ESI HD]. The recent cosmological 
observations from Planck2013 impose tight bounds on the EoS parameters in an emergent 
universe (EU) |37j . For a viable cosmological scenario, it is further necessary to consider 
a consistent model of the universe which contains radiation dominated, matter dominated 
and subsequently the late accelerated phases of the universe. In the original EU model 
it is shown |2nj that the composition of the universe is hxed once A is hxed. A problem 
thus arises as to how a pressureless matter component could be accommodated within such 
a scenario. However, allowing interaction among the constituent fluids of the emergent 
universe may open up richer physical consequences. 

Interacting fluid cosmological models have been previously considered in the literature, 
and among a variety of reasons and motivations for such models, analyses of interactions 
within the dark sector are quite popular |41l 1421145114411451146] . In the present context in¬ 
teraction among the constituent fluids is useful to obtain a consistent evolutionary scenario 
of the universe. Another important consistency condition is imposed through the thermo¬ 
dynamics of an expanding universe. There has been a lot of recent interest in the study 
of the connection between thermodynamics and gravitational dynamics in the presence of 
horizons [izlIlslIlgllsnilSIlliaiMl El ES] . Another aspect the Hawking temperature of the 
apparent horizon in a FRW universe was determined by Cai et al. |56j . Indeed, Einstein’s 
equations have been interpreted as a thermodynamical relation resulting from the displace¬ 
ment of the horizon. Since the emergent universe scenario entails a phase of accelerated 
expansion, it is relevant here to study the status of the second law of thermodynamics in 
the picture involving interacting fluids in the emergent universe. 

With the above motivations, in the present paper we consider two different cases in the 
emergent universe scenario: (i) a two-fluid model with interaction of the fluid having the 
non-linear EoS given by eq. (1) with another barotropic fluid, beginning at some time t = ti 
(Model-I) and (ii) a three-fluid model with interaction among the various constituent fluids 
with different individual equations of state, starting at a time t = to (Model-H). The paper 
is organized as follows : In Sec. 2 we set up the field equation for emergent universe in 
the general theory of relativity, and obtain cosmological solutions. In Sec. 3 we consider 
the above two cases, Model-I and Model-II. We show that generation of pressureless matter 
fluids is possible within the emergent universe scenario. In Sec. 4, we demonstrate the 
consistency of the emergent universe scenario with interacting fluids with the generalized 
second law of thermodynamics. Finally, we make some concluding remarks in Sec. 5. 

2 Field Equation and Cosmological Solution 

The Einstein held equation is given by 


where Rf^u, R, 9fiiy, T^u and G represent the Ricci tensor, Ricci scalar, metric tensor, matter- 
energy tensor and Newton’s gravitational constant. Here we consider four dimensions for 
which /i, v runs from (0, 1, 2, 3). 

We consider a flat Robertson-Walker metric which is given by 


= —di^ + a^{t) dr'^ -|- + sin^OdcfP'^ 


(3) 


where a{t) represents the scale factor of the Universe. We consider the energy momentum 
tensor as Tjf = diagonal{p, —p, —p, —p), where p is the energy density and p is the pressure. 
Using the flat Robertson-Walker metric given by eq. (3) in the Einstein’s field equation one 
obtains 

(4) 

(5) 




p= 


a 

2 - + 
a 


where we set G = -^ and c = 1. The conservation equation is given by 


dp 

dt 


+ ^H {p + p) = Q 


( 6 ) 


where H = ^ represents the Hubble parameter. As mentioned earlier, Mukherjee et al. 
|20j obtained an emergent universe scenario with a polytropic equation of state (henceforth, 
EoS) given by 

p = Ap — Bp^/"^ (7) 

where A and B are arbitrary constants. Making use of the conservation equation and the 
EoS given by eqs. (4) and (5) in eq. (7), one obtains a second order differential equation 
given by 

2--7 (3A-7 1) -\/3R- = 0. (8) 

a \aj a 

The scale factor of the universe is thus obtained integrating eq. (8) which is given by 


a{t) 


'3K{A + l) 
2 




2 

3(A+1) 


(9) 


where K and a are the two integration constants. It is interesting to note that B < 0 
leads to a contracting universe whereas with B > 0 and A > — 1 leads to a non-singular 
solution which is expanding. The later solution corresponds to an emergent universe which 
was obtained by Mukherjee et. al. [20]. The energy density of the universe in terms of scale 
factor is obtained from eq. (6) making use of EoS given by eq. (7) which is given by 


p{a) 


1 

(A+ 1)2 



K 


3(A+1) 

a 2 


2 
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Expanding the above expression, one obtains energy density as the sum of three terms 
which can be identified with three different types of fluids. Thus, the components of energy 
density and pressure can be expressed as follows: 


p{a) = T^i^iPi and p{a) = 


( 11 ) 


A 














where we denote 


52 2 KB 1 K2 1 

“ (^+ 1 ) 2 ’ ~ (yl+l)2^3(^> P3 - + 1)2 q3(A+1) ( 12 ) 

52 ^ - 1) 1 _ AK"^ 1 

"(yl + l) 2 ’ “ (yl + 1)2 “ (A + 1)2 a3(^+i)' 

Comparing with the barotropic EoS given by pi = uiiPi one obtains cji = — 1, a ;2 = and 
UJ3 = A. Thus, the parameter A plays an important role in determining the composition of 
the fluids in the universe. For example, A = ^ leads to a universe with radiation, exotic 
matter and dark energy, A = Q leads to dark energy, exotic matter and dust. Thus once 
the EoS parameter A is fixed the composition of the fluid in the universe gets determined. 
In order to obtain a viable scenario of the universe we consider interacting fluids model in 
the next section so that a transformation of one kind fluid in the later epoch gives rise to 
the composition of matter that we observe today. 


3 Cosmological Models 

In this section we consider two different models of interacting fluids in an EU scenario. 

Model I : The two fluids model 

In this case, we consider two interacting fluids with densities p and p' respectively which 
can exchange energies with each other. One of the fluid with energy density, say p is 
dominated to begin with satisfying a non-linear EoS given by eq. (1) which leads to an 
emergent universe model as discussed above with no interaction. The contribution of the 
other fluid in the energy density of the universe is assumed to be important at a later epoch. 
The corresponding pressure of the former fluid is given by 

p = Ap — Bp^!"^ . (14) 

where A and B are constants. The other fluid satisfies a barotropic equation of state which 
is given by 

p' = (jj'p' (15) 

where oj' corresponds to EoS parameter. The Hubble parameter (^H = ^ evolves according 
to the Friedmann equation which is given by 

3 Er 2 = p + p'. ( 16 ) 

In this case we consider a cosmological model where exchange of energy between two different 
fluids is allowed.The are many astrophysical and cosmological motivations for considering 
energy exchanges between the various components of the universe [39-44]. Different phe¬ 
nomenological considerations dictate the onset of such interactions. For example, in various 
scalar field models of dark energy, such as in quintessence, or k-essence, there arise phase 
transitions during particular eras resulting in decay of the cosmological vacuum energy, as 
well as particle production. Similarly, there could be other cases of energy exchange, such 
as due to the evolution of a population of primordial black holes whose evaporation time 
depends on the particular formation mechanism or formation era. In the present paper. 



without assuming any specific mechanism for energy exchange, we assume that the inter¬ 
action starts at some particular time ti. The two interacting fluids respect a total energy 
conservation equation and their densities evolve with time as 


p + ‘iH{p + p) = -apH, 

(17) 

p' + 3Hip' + p’) = apH 

(18) 


where a represents a coupling parametrizing the energy exchange between the fluids. One 
may view the above interaction as a flow of energy from first kind of fluid to the second one 
(say, dark matter) beginning at the epoch considered here. Now, using eq. (14) in eq. (17), 
we get a first order differential equation which can be integrated to obtain the behaviour 
of energy density in terms of the scale factor of the universe and the interaction coupling 
factor. Thus the energy density and pressure for the fluid of the first kind are given by 
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{A + l + ^y 
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2KB 


(^ + 1+^)2 3(A+l+j 
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(A + 1 + f )2 a3(^+i+f) 
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KB{A-1 + ^) 1 
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(7l + l + f)2 (24+1 + «)2 3(A+l+f) (^+1+ a)2^3(A+l+|) 


( 20 ) 


If the interaction is with a pressureless dark fluid i.e., p' = 0 (however, p' / 0), eq. (18) 
can be integrated using eqs. (16) and (19) which determine the total energy density and 
pressure as follows : 


, _ ^2 2KB 1 772 1 

Ptotal -P + P - + A + 1 a3(A+l) ’ 

^2 KB{A-1) 1 ^ 24772 I 

Ptotal -P- + (^ + 1)2 + 1)2 ^3(A+1) ' 

The equation of state parameter for the second fluid is given by 

/ P Ptotal P 

U1 = - = -. 

P ptotal P 


( 21 ) 

( 22 ) 


(23) 


In the limiting case as u' —0, we get ptotai = P- An interesting case emerges when the 
coupling parameter a = 2 and A = ^. In this case a universe with dark energy, exotic 
matter and radiation to begin with (i.e., before the interaction sets in) made a transition 
to a matter dominated phase. Thus, an emergent universe with radiation, dark energy 
(component that is subdominant at early times) and exotic matter (for A = |) to begin 
with transits to a universe with matter domination phase after an epoch after t > ti for the 
interaction coupling strength a = 2. Hence, a consistent scenario of the observed universe 
in the EU model may be realized in this case. 


Model II : The three fluids model 


The original EU model was obtained in the presence of non-interacting fluids permitted 
by the parameter H in a flat universe case. The corresponding densities and pressures are 
given by eqs. m and (fT3]l respectively. Eor non-interacting fluids, the EoS parameters 
for the three fluids permitted above are given by ioi = —1, UI 2 = \ {A — 1), = A. Eor 

0 < H < 1, it accommodates dark energy, exotic matter and the usual barotropic fluid. The 
energy density and pressure of the exotic matter and that of the barotropic fluids decreases 



















with the expansion of the universe. However, the rate of decrease is different evident from 
eqs. (m and ( fT3]) . We assume an interaction among the components of the fluid in the 
universe which is assumed to be originated at a later epoch (such interactions could arise 
due to a variety of mechanisms [m sa SHI El SSI SB]). Assuming onset of interaction 
among the composition of the fluid at t > to, the conservation equations for the energy 
densities of the fluids now can be written as 


Pi + 3H{pi + pi) — —Q', (24) 

P 2 + 3H{p2+P2) = Q, (25) 

P3 + ^H(p3 + Ps) = Q' — Q, (26) 

where Q and Q' represent the interaction terms, which can have arbitrary form, pi represents 
dark energy density, p 2 represents exotic matter and ps represents normal matter. In 
this case Q < 0 corresponds to energy transfer from exotic matter sector to two other 
constituents, Q' > 0 corresponds to energy transfer from dark energy sector to the other 
two fluids, and Q' < Q corresponds to energy loss for the normal matter sector. The case 
Q = Q' corresponds to the limiting case where dark energy interacts only with the exotic 
matter. It is important to see that although the three equations are different the total 
energy of the fluid satisfies the conservation equation together. It is possible to construct 
the equivalent effective uncoupled model, described by the following conservation equations: 


Pi + 3H{1 + ool-^-^)pi = 0 

(27) 

P 2 + 3H{1 + 002 ^ ^)p 2 = 0 

(28) 

P 3 + 3H{1 + ioJ^)p 3 = 0 

(29) 


where the effective equation of state parameters are given below: 


UJ 


eff _ 


= 001 + 


Q' 

3Hpi' 


(30) 


eff 

OJ 2 =002 — 

eff I 

OO^-' •' = W 3 + 

Now, if we consider the interaction as Q — Q' = 
the normal fluid becomes 


Q' 

3Hp2' 

Q-Q' 

3Hp3 

= —/3Hp^, the effective state parameter for 






eff _ 


= OJ 3 


3 


(33) 


In fig. 1 we plot the variation of effective equation of state parameter 00 ^^^ with 003 (which 
corresponds to A of the EoS parameter) for different strengths of interaction determined 
by /?. We note that as the strength of interaction is increased the value of 0:3 (i.e., A) for 
which oo^ff = 0 (corresponds to matter domination) is found to increase. Thus a universe 
with any A value is found to admit a matter dominated phase at a late epoch depending 
on the strength of the interaction which was not permitted in the absence of interaction 
in an EU model proposed by Mukherjee et al. [2^. It may be pointed out here that in 
the very early era a universe is assumed with a composition of three different fluids having 
no interaction in this picture, thus the behaviour of the universe at early times remains 
unchanged as was found in the original EU model. Thus the emergent universe scenario 
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Figure 1: Plot of with EoS parameter A for different interaction (3. The thick, dash, 
red line and green lines are for (3 = 0.5, 1, 2, 3 respectively. 


proposed by Mukherjee et al. m can be realized in the early era but at a later epoch the 
composition of matter changes in the present scenario from its original one with the onset 
of interaction. This feature represents a clear improvement over the earlier cosmological 
scenario in an emergent universe [2nii23En] where it is rather difficult to accommodate a 
pressureless fluid. 

4 Generalized Second Law of Thermodynamics in Emergent 
Universe Model 

In the above analysis emergent universe models have been discussed in the presence of differ¬ 
ent types of interaction among the fluids. In the present section we investigate consistency 
of thermodynamic properties considering the universe as a thermodynamical system. For a 
flat geometry the apparent horizon coincides with the Hubble horizon which is given by 

^ (34) 

In general, the apparent horizon is a function of time. Thus, a change in the apparent 
horizon leads to a change in volume, and consequently, the energy and entropy will change 
by dE and dS respectively. The energy momentum tensor before and after the change are 
described by the same T^,y, and we can consider that the pressure and the temperature 
remains the same HU sg ig [501 El]. The first law of thermodynamics for the fluids 
considered here are given by 

dSi = ^ {pidV + dEi) 
dS2 = ^ {p2dV + dE2) 
dSs = ^ {p3dV + dEs) 




( 35 ) 





where the volume of the system V = is bounded by the apparent horizon, and thus 
dV = ATTr^dr^. Hence, the rate of change of entropy for the above fluids are given by 

111 
<5i = ^ {A'n:r\fAVi + ^i) , *^2 = ^ {A.'Kr\rAP 2 + -^ 2 ) , Ss = — (Anr^fAPs + ^ 3 ) (36) 

where 

= 2^a ((1 + ^i)pi + (1 + ^2)P2 + (1 + <^ 3 )^ 3 ) (37) 

which is obtained differentiating the equation 

— = \{pi+P2 + ps) (38) 

'A 'J 


and then making use of the conservation equations given by eqs. (27)-(29) with eq. (34). 
The corresponding energy density and pressure are given by 

^ dvr o 

^1 = 

^ dvr o 

^2 = Y^>2, 

47r Q 

^3 = (39) 

and 

eff eff eff 

Pi=a;/Vi, P2=(^2 P2, P3 = ^^3 P3- (40) 

Using the time derivative of the apparent horizon we get 

51 = - HrA){l + Wi^^) 

5 2 = ^^^P2{rA - HrA)il + 

53 = ^^P3{rA- HrA)il + uj^^^). (41) 

According to the generalization of black hole thermodynamics [501ICT [52] to a cosmological 
framework, the temperature of the horizon is related to its radius [5311541 ?] as 


Th 


1 

27rrA ’ 


leading to the rate of change of entropy given by 


Sh = lOvr^r^rk. 


(42) 


(43) 


It may be mentioned here that the Hawking radiation of apparent horizon in a FRW universe 
is computed by Cai et al. [56] . The rate of change of total entropy becomes Stotai = 
Si + S 2 + S 3 + Sh, which leads to 

Stotai = A-K‘^r\H [(1 + a;i)pi + (1 + i02)p2 + (1 + ^^ 3 )^ 3 ]^ > 0. (44) 

The non-negativity of the time rate of change of Stotai demonstrates the validity of the 
second law of thermodynamics in the context of the emergent universe model. 


q 







5 Discussion 


In this paper we have investigated cosmology of the emergent universe scenario in the 
presence of interacting fluids. The purpose of the present analysis is to demonstrate the 
possibility of obtaining viable cosmological dynamics of the emergent universe. Two dif¬ 
ferent cosmological models have been presented here. In Model I, we consider the flow of 
energy from the fluids required to realize the emergent universe to a pressureless fluid which 
sets in at an epoch t = ti. The density of the pressureless fluid assumes importance as mat¬ 
ter component after the epoch ti. In Model II, we consider interactions among the three 
fluids of the emergent universe at time t = to. Before this epoch the emergent universe can 
be realized without an interaction among the fluids. The problem with earlier cosmological 
realizations of the emergent universe was that once the EoS parameter A is fixed at a given 
value, the universe is unable to come out of the phase with a given composition of fluids. In 
the present work we overcome this problem by assigning an interaction among the fluids at 
the epoch to. A cosmological evolution of the observed universe through unified dynamics 
of associated matter and dark energy components thus becomes feasible in the emergent 
universe scenario. In fig. (1), we plot variation of with EoS parameter cus = A for 
different interaction. It is evident that a early universe with a radiation dominated phase 
transits to a matter dominated phase with all the features observable at the present moment 
with the onset of interaction considered here. Consistency of the interacting fluid emergent 
universe scenario with the generalized second law of thermodynamics is also shown. Eur- 
ther work is needed in order for a comparative analysis of the interacting fluid emergent 
universe cosmology with more popular current cosmological models. Detailed analysis of 
observational constraints pertaining to various eras of the universe are expected to yield 
bounds on the parameters of the emergent universe models considered here, leading to a 
firmer assessment of the viability of such models. 

6 Acknowledgement: 

BCP would like to thank S. N. Bose National Centre for Basic Sciences, Kolkata for warm 
hospitality during a visit when the work was initiated. BCP likes to thank Institute of 
Mathematical Sciences (IMSc), Chennai for hospitality during a visit. BCP also acknowl¬ 
edges University Grants Commission (UGC), New Delhi for financial support awarding with 
a Major Research Project (Grant No. P. No. 42-783/2013(SR), dated 13 Feb 2014). 



References 

[1] Penzias A A and Wilson R W 1965 Astrophys. J. 142 419 

[2] Guth A H 1981 Phys. Rev. D 23 347 

[3] Sato K 1981 Mon. Not. Roy. Astron. Soc. 195 467 

[4] Linde A D 1982 Phys. Lett. B 108 389 

[5] Albrecht A and Steinhardt P 1982 Phys. Rev. Lett. 48 1220 

[6] Albrecht A 2000 arXiv:astro-ph/0007247 

[7] Linde A D 1983 Phys. Lett B 129 177 

[8] Linde A D 1985 Phys. Lett B 162 281 

[9] Mukhanov V F and Chibisov G V 1981 JETP Lett 33 532 

[10] Riess A G et al. 1998 Astron J. 116 1009 

[11] Perlmutter S et al. 1998 Nature 51 391 

[12] Perlmutter S et al. 1999 Astrophys. J., 517 565 

[13] Garroll S M 2001 Living Rev. Rel. 4 1 

[14] Dicke R H, Peebles P J E, Roll P J and Wilkinson D T 1965 Astrophys. J. Lett. 142 
414 

[15] Ellis G F R, Maartens R 2004 Class. Quant. Grav. 21 223 

[16] Harrison E R 1967 Mont. Not. Roy. Aston. Soc. 69 137 

[17] Ellis G F R, Murugan J and Tsagas C G 2004 Class. Quant. Grav. 21 233 

[18] Murlryne D J, Tavakol R, Lidsey J E and Ellis G F R 2005 Phys. Rev. D 71 123512 

[19] Padmanabhan T and Padmanabhan H 2014 Int. J. Mod. Phys. D 23, 1430011 

[20] Mukherjee S, Paul B C , Dadhich N K, Maharaj S D and Beesham A 2006 Class. 
Quant. Grav. 23 6927 

[21] Ghimento L 2004 Phys. Rev. D 69 123517 

[22] Mukherjee S, Paul B C, Maharaj S D and Beesham A 2005 arXive:gr-qc/0505103 

[23] Paul B C and Ghose S 2010 Gen. Rel. Grav. 42 795 

[24] Banerjee A, Bandyopadhyay T and Ghakraborty S 2008 Gen.Rel. Grav. 40 1603 

[25] Banerjee A, Bandyopadhyay T and Ghakraborty S 2007 Grav. & Cosmology 13 290 

[26] Debnath U 2008 Class. Quant. Grav. 25 205019 

[27] del Gampo S, Herrera R, Labrana P 2007 J. Cosmo. Astropart. Physics 30 0711 

[28] Beesham A, Chervon S V and Maharaj S D 2009 Class. Quantum Grav. 26 075017 



[29] Chervon S V, Maharaj S D, Beesham A and Kubasov A 2014 An emergent universe sup¬ 
ported by chiral cosmological fields in Einstein Gauss-Bonnet Gravity arXiv: 1405.7219 

[30] Beesham A, Chervon S V, Maharaj S D and Kubasov A 2013 Quantum Matter 2 388 

[31] Beesham A, Chervon S V, Maharaj S D and Kubasov A 2009 Glass. Quantum Grav. 
26, 075017. 

[32] Bag S, Sahni V, Shtanov Y and Unnikrishnan S 2014 JGAP 07 034 

[33] Majumdar A S 2001 Phys. Rev. D 64 083503 

[34] Bose N and Majumdar A S 2009 Phys. Rev. D 79 103517 

[35] Bose N and Majumdar A S 2009 Phys. Rev. D 80 103508 

[36] Komatsu E et al. 2011 Astrophys. J. SuppL, 192 18 

[37] Marra V, Amendola L, Sawicki I, Valkenburg W 2013 Phys. Rev. Lett. 110 241305 

[38] Paul B C, Thakur P and Chose S 2010 Mon. Not. Roy. Astron. Soc. 407 415 

[39] Paul B C, Chose S and Thakur P 2011 Mon. Not. Roy. Astron. Soc. 413 686 

[40] Chose S, Thakur P and Paul B C 2012 Mon. Not. Roy. Astron. Soc. 421 20 

[41] Barrow J D and Clifton T 2006 Phys. Rev. D 73 103520 

[42] Chimento L P 2010 Phys. Rev. D 81 043525 

[43] Jamil M, Saridakis E N and Setare M R 2010 Phys. Rev. D 81 023007 

[44] Lip S Z W 2011 Phys. Rev. D 83 023528 

[45] Costa F E M, Alcaniz J S, Jain D 2012 Phys. Rev. D 85 107302 

[46] Cotsakis S and Kittou C 2013 Phys. Rev. D 88 083514 

[47] Frolov A V and Kofman L 2003 J. Gosmol. Astropart. Phys. 05 009 

[48] Danielsson U H 2005 Phys. Rev. D 71 023516 

[49] Bousso R 2005 Phys. Rev. D 71 064024 

[50] Cai R C and Kim S P 2005 JHEP 02 050 

[51] Akbar M and Cai R C 2007 Phys. Rev. D 75 084003 

[52] Cibbons C W and Hawking S W 1977 Phys. Rev. D 15 2738 

[53] Jacobson T 1995 Phys. Rev. Lett. 75 1260 

[54] Padmanabhan T 2005 Phys. Rep. 406 49 

[55] Paranjape A, Sarkar S, Padmanabhan T 2006 Phys. Rev. D 74 104015 

[56] Cai R C, Cao L M and Hu Y P 2009 Glass. Quantum Grav. 26 155018, arXiv;0809.1554 


